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We introduce a deductive statistical mechanics approach for granular materials which is for- 
mally built from few realistic physical assumptions. The main finding is an universal behavior for 
the distribution of the density fluctuations. Such a distribution is the equivalent of the Maxwell- 
Boltzmann's distribution in the kinetic theory of gasses. The comparison with a very extensive set 
of experimental and simulation data for packings of monosized spherical grains, reveals a remark- 
ably good quantitative agreement with the theoretical predictions for the density fluctuations both 
at the grain level and at the global system level. Such agreement is robust over a broad range of 
packing fractions and it is observed in several distinct systems prepared by using different methods. 
^ , The equilibrium distributions are characterized by only one parameter (fc) which is a quantity very 

O ' sensitive to changes in the structural organization. The thermodynamical equivalent of k and its 

' relation with the 'granular temperature' are also discussed. 
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I. INTRODUCTION 



^^-H , Granular materials are central in a very wide range of domains from agriculture to pharmaceutical industry, but, 
^ ' despite such a central role in most fields of human activity and their ubiquitous presence in scientific research areas, 
^ ; _ their behaviour and properties remain elusive. Granular materials differ from all the other materials: they can flow 
. like liquids under some conditions and act as solids under others [J, [2j @] ■ Their intrinsic complexity is consequence of 
^ ■ the fact that they are composed of many pieces (the grains) that can assemble in large structures behaving collectively. 
I , The study of granular materials has forced scientists to rethink established classifications of matter and to reformulate 
' O statistical mechanics in a new context. These studies are springing out new ideas that are advancing the understanding 
S , of large classes of complex materials from composites to foams. 

Q ■ The question whether a statistical mechanics approach can correctly describe granular swtems has sparkled an 
lZIi ' intense debate in the scientific community in recent years 0, d, i, 0, B i, [3, [H [ll H Q d, [la, [13, d, E S [SI, ill . 
Such a statistical mechanics theory of granular matter must be able to connect the observable behavior of large 
aggregates with the structural and dynamical properties of the grains constituting these systems. The great interest 
in such 'thermodynamical' theory of granular packings resides in the fact that it would lead to a description of such 
a complex system in terms of a few relevant parameters only 0, [1, [^, [isl . [isj . This is in analogy with traditional 
, thermodynamics, where the state of a large system comprised of many interacting elements can be completely described 
' by knowing temperature, volume or pressure. 
' \ In this paper we show that such approach is formally possible: a statistical mechanics formulation for granular matter 
I ■ can be built deductively from a small set of physical assumptions and the equivalent of the Maxwell-Boltzmann's 
f — ] distribution can be formally derived. More significantly, we show that such a distribution describes very accurately 
. empirical and numerical data concerning volume distributions in sphere packings prepared with several different 
^ methods in a broad range of packing fractions. 
• i-H , One of the characteristics that differentiate granular materials from molecular materials is that frictional forces 
^ ■ among grains can dissipate energy and drive the system towards frozen, jammed configurations. Because of such 
^ \ energy dissipation, granular materials can change their static configurations only when energy is injected in the 
■ - - 1 system. This can be done in several ways, for instance by tapping the system with vertical vibrations or by pulsing 
fluids flows in fluidised beds or by shearing or by rotating the container. The relevant fact is that such action should 
make possible for the system to change configuration and explore the accessible phase-space. 

The statistical mechanics approach proposed in this paper is built within the same framework of the deductive 
foundation of statistical mechanics proposed by O. Penrose in [2^. The present paper is structured as follows. In 
Section ini we introduce the general idea of the system's observational states and the probability of transition between 
states when the system is driven by trains of energy injections. In Section [llll we define a microscopic 'fine grained' 
description of system's state concerning grain packs at static equilibrium. In this section we also discuss other possible 
'coarse grained' macroscopic descriptions and the links between these different descriptions. A classification of states 
in terms of their transitional probabilities is proposed in Section HVl where the concept of ergodic set is also introduced. 
In Section |Vl we discuss the statistical equilibrium. The system entropy is defined in Section IVTl where a maximum- 
entropy prediction for the probability of a given (coarse-grained) state is also introduced. The readers that are only 
interested on the theoretical outcomes might skip all the previous sections and jump directly to Section IVIII where 
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the probability density function for the system's volume fluctuations is calculated. In Section IVIIIl the concept of 
'granular temperature' and the links with volume fluctuations at equilibrium are discussed. In Section HXl it is shown 
that the theoretical predictions for the probability density function for local volumes at the grain levels (Voronoi' 
cells) reproduce extremely well experimental and simulation data. Moreover, we show that the theory predicts also 
the global fluctuations on the scale of the whole system. The signatures of transitions occurring at the random 
loose packing limit and at the random close packing limit are highlighted by studying the changes occurring at these 
densities in a quantity (fc) which is shown to be the analogous to the specific heat for these systems. Conclusions and 
perspectives are given in Section [Xj 



II. OBSERVATIONAL PROBABILITY AND TRANSITION BETWEEN STATES 



Let us here consider a granular system that evolves by means of repeated energy injections. We are interested in 
the system's state observed in between trains of energy injections when all the grains are at rest in a stable static 
configuration. Let us use the term 'trials' to indicate such energy injection trains. Let us define Pt(^' ^ ^'o) as the 
probability to find the system in the observational state \1/ after an experiment consisting of t trials, given that at 
the beginning of the experiment {t = 0) the system was prepared in the state ^q. We will return on the definition of 
the observational state in the next section. In this section we assume that the system static state can be defined and 
encoded. 

A further important assumption is that the system is Markovian (23j : the current observational state embodies 
all the observational information about the past history of the system that are relevant to its future observational 
state. Such assumption allows us to write the Chapman-Kolmogorov equation describing the evolution under energy 
injections 23]: Pt+i('^ ^ ^'o) = ^ ^i)-ft(^i ^ ^o); which can be re- written in the more familiar form 

= ^u;(*^^'i)p(^'i,i) ; (1) 

where ^ = Pi(^' ^ 'i'l) is the probability of transition between the state "ifi to the state by performing 
one trial on the system. Whereas, p{'^i,t) is a simplified notation for Pt(^'i, 4'o) with the initial state omitted. 

The Markovian assumption is quite strong from a statistical point of view. However, in the context we are in- 
vestigating, this assuption is adequate. Indeed, a granular system at rest in a static stable configuration can stay 
'frozen' in such a state for ever. Such a state might depend on all the system past history, the Markovian hypothesis 
simply assumes that all the relevant information which will determine the system future evolution is contained in the 
properties of the present system's state. Indeed, there are no simple mechanisms to store 'hidden' information about 
the granular system which are not embodied into the complete microscopic description of the system at rest. It must 
be noted that, on the contrary, a coarse-grained description of the system's state might classify as identical states 
with different microscopic properties. Although belonging to the same coarse-grained classification, such states are 
distinct at microscopic level and might evolve with different transition probabilities. 



III. DEFINITION OF THE SYSTEM'S STATE 



The aim of this section is to show with a practical example that it is formally possible to precisely encode the system's 
microscopic state. Granular materials are made of several grains which are typically confined in a container. Each 
grain 'i' in the system can have a different shape (s^) and a different mass (m^). The geometrical disposition of each 
grain can be characterized by the coordinates of its barycenter = {xi, yi, Zi) and the grain orientation = (9i, 4>i). 
If we assume that: 1) the grains are non-deformable objects but they can overlap; 2) the grains interact only when 
they are in contact or they are overlapping; 3) the interaction energy is a known function of the grain overlaps; then, 
it is possible to give a complete description of the observable state of a packing of N grains at rest in a stable static 
configuration by providing the set of positions and orientations of all grains: ^' = {ri, r2, r^, Oi, 02, o^r}. 

However, when friction is present, the tangential frictional force between two grains is a function of the relative 



path that has produced such static configuration [2j|. Such dependence on the previous history would falsify the 



Markovian assumption if we adopt the previous classification for 5*. On the other hand, we can solve this apparent 
difficulty by including the set of tangential forces between grains into the microscopic description. In this way all the 
information about the relevant past history is encoded in the preset state description. The future evolution depends 
only on the present state and the trials. The theory satisfies therefore the Markovian assumption. Such an encoding 
is for instance used in numerical simulations of granular systems such as the Discrete Element Method (DEM) with 
'soft' spheres [25, 26, 21]. 



3 



Independently on the details of the system and on its specific encoding, the previous example has shown a very 
important fact: the system microscopic state can be completely encoded in terms of a finite set of quantities. However, 
we must consider that there are several possible descriptions of the system's state and the use of different variables 
might be more appropriate in different systems. For instance, it can be convenient to include in the local description 
of the packing also some information about the surrounding configuration. To this purpose one can subdivide the 
packing into a set of k local 'elementary cells'. Such cells can be viewed as the portions of space associated to local 
configurations in the granular pack. If we call c the properties of such cells (size, shape, volume, connectivity, center of 
mass, orientation, etc.), then the system's state can be conveniently classified in terms of the properties of such local 
'elementary cells': ^E* — {ci, C2, Cfc}. Let us note that the two descriptions in terms of positions and orientations 
or in terms of local cells must be equivalent and interchangeable. The complete information about the system's 
configuration is encoded in both the representations. In the following we will refer to the state 4" as a 'microscopic' 
state. 

The classification in terms of the microscopic state ^ is the most complete description of the system's static 
configuration. However, in most of the practical cases such a 'fine grained' description is neither necessary nor useful. 
There are indeed several different microscopic configurations that share the same macroscopic physical properties. 
A 'coarse grained' description of the system's state can be therefore more useful, and practical, when one wants to 
associate structural properties to physical properties. An example of a macroscopic - coarse grained ~ description of 
the system's state is its classification in terms of the total volume V occupied by the system (or equivalently in terms 
of the packing fraction, volume of the grains / V). Certainly, there is a very large number of microscopic states ^' 
which correspond to the same volume V. Nonetheless, there are several structural and physical properties of these 
systems that can be unambiguously classified in terms of the packing fraction only [2^ . The possibility to construct 
a coherent framework which allows one to consistently describe the system at different levels of detail is one of the 
main features of statistical mechanics (23j . A statistical mechanics theory of granular material should provide the 
instruments to unambiguously describe the average system properties and their probability distributions at any chosen 
level of description. 



The system's states can be classified in terms of their transitional probabilities f2J||. In particular, a state $ is 
transient if there exists a state 5* that the system can reach from $ but it cannot go back to $. If the state is 
not transient it is said to be persistent. The set 2'($) of all states that can be reached from a persistent state <& 
is called ergodic set. Persistent states may be separated into one or more distinct ergodic sets. Once the system 
reaches a persistent state it will remain in the ergodic set containing such state and it will repeatedly visit all the 
states in the set never going outside. For large t the system will eventually reach the limit probability distribution 
p(^', t 00) ~ poo(*), that is the solution of the system of equations Pooi^) — X)* "^(^ ^ ^)j'oo(^)- Note that such 
limiting probability is stationary in the sense that it is not dependent on t. A system which is confined in an ergodic 
set Z and visits the states € Z with frequencies proportional to Poo(^') is said to be at statistical equilibrium. 

We note that such transitional probabilities, which define the classification of states and ultimately the system's 
statistical equilibrium, are intrinsically dependent on the kind of trials we are performing on the system. Therefore, 
the equilibrium state and even the ergodic set are not solely associated with the kind of system in exam. Conversely, 
they are uniquely associated to the combination of both system properties (grain shapes, roughness, friction, container 
shape, ... ) and the kind of trial (tapping, fluid flow, pouring, rotating, ...) and its properties (vibration frequency 
and amplitude, fluid intensity, pouring height, rotation velocity,...). 

These classifications of the system's states and all the transition properties and evolution relations (Eq. [T|) discussed 
in this section and in the previous one, can be applied also to coarse grained states, providing that they satisfy the 
Markovian assumption. It is of some interest to note that a coarse grained state A can be at the statistical equilibrium 
even if the microscopic states haven't reached equilibrium yet. 



We have already introduced in the previous section the definition of statistical equilibrium. This is a very important 
concept that must be discussed with some care. In a molecular gas, statistical equilibrium is reached after a certain 
time when the system has achieved a homogeneous internal temperature and it has thermalized to the surrounding 
environment. In the context of granular materials and in particular in tapping experiments, the 'equilibrium' state 
is typically considered reached when, after a certain number of tappings with a given intensity, the packing fraction 
begins to oscillate around a constant value that is independent from the starting state. Such an heuristic definition is 



IV. CLASSIFICATION OF SYSTEM'S STATES 




V. STATISTICAL EQUILIBRIUM 
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consistent with the formal one given in the previous section, however there can be cases when more than one ergodic 
set can be reached from the same starting configuration. Each one of such ergodic sets can have different statistical 
equilibrium properties. 

It must be noted that the definition of statistical equilibrium is intrinsically associated with the transitional prob- 
ability and therefore with the kind of trial that we are repeatedly performing on the system, it is not a property of 
the system's static configuration. In this sense it is misleading to think that statistical equilibrium can be reached 
by means of certain prescribed operations only (such as tapping or fluid flows). On the contrary there are a very 
large number of actions (trials) that will drive the system toward equilibrium. For instance, if we pour grains from a 
vessel into another and we do it several times being careful to repeat the operation each time in the same way (same 
height, inclination angle,...) we will eventually reach a situation in which the static packing configurations have stable 
statistical properties. Analogously, statistical equilibrium can be achieved by rotating the container or by blowing 
air from below or by stirring the aggregate or by several other possible manipulations. The open question is whether 
these 'equilibrium' states are the same states or not. They might have the same average packing fraction but, to 
coincide they must also belong to the same ergodic set. Generally speaking, it is quite unlikely that exactly the same 
regions of the phase-space are explored by means of operations as different as fluidized beds or pouring. However, 
we will show in Section IIXI that there are evidences of remarkable universal properties that hold across very diverse 
preparation methods and different numerical simulations. Such universality is suggesting that there might be a very 
large class of processes that leads to the same ergodic set or to sets that are mostly overlapping. 

VI. MAXIMUM ENTROPY PREDICTION FOR THE STATE PROBABILITIES 

There are several different definitions of Entropy that can be applied to different contexts such as thermodynamics 
of molecular gasses or information theory. In this paper we will use the 'statistical entropy' (or Gibbs entropy) S{Z) 
which is better suited to be applied to the study of the ensemble of states Z constituting the ergodic set: 

S{Z)^-Y,Poo{A)\ogpoo{A)+Y,Poo{A)S{A) . (2) 
A£Z Aez 

Here the state A is a given 'coarse grained' description of the system; Poo{-A) is its - stationary - probability within 
the ergodic set Z and S{A) is the entropy of the state A, which is: S{A) = — J2<sreA '?(^) log where <7(\E') are the 
a-priori probabilities to find the system in the microscopic state £ A. When all the microscopic states are a-priori 
equi-probable, the expression for S{A) becomes 

S{A)^logn{A) , (3) 

with Q,{A) the number of microscopic states which are classifiable under the same coarse grained state A. Eqsl2]and 
[3] give a coherent definition of the system's entropy independently on the level of details at which the system's states 
are described. In general, the choice for the tuning of the system's description depends on our experimental capability 
of measuring the structural properties, and on our analytical or computational ability to evaluate the number of 
microscopic configurations which are contributing to such states. For what concerns the a-priori probabilities (?(^'), 
in the following we will assume that, a-priori, the probabilities to find the system in any given microscopic state are 
all equal and therefore we will use Eq[3) This is not an essential assumption for the present theory, however this is 
the most likely scenario for a granular system prepared by means of some random mixing operation. 

Let us adopt characterize the system's state in terms only of its total occupied volume V (i.e. A = V). Different 
experiments can result in granular packings with different total occupied volume V. For instance, in the classical 
experiment by ^29i . |3Q . iSli] . different average packing fractions are obtained by driving the system with different 
tapping intensities. Here we are interested in the properties of the ergodic set, which is the set of all the possible total 
volumes which can be reached by means of the chosen system's handling. The probability of each of such realizations 
is Pociy)- The ergodic set at statistical equilibrium is characterized by the average volume V (or equivalently the 
average packing fraction p) and the study of the statistical properties of the ergodic set coincides with the study of 
the volume fluctuations around such average values. 

The most important property of entropy is that it is a non decreasing function, therefore at the statistical equilibrium 
the entropy must reach its maximum. We can calculate PociV) by searching for the functional form of the probability 
distribution function which maximizes the entropy S{Z) (EqU). Such maximization must be done under the condition 
that the average occupied volume is equal to V, which is the characteristic average volume associated to the ergodic 
set Z corresponding to the specific trial. This yields to: 

Poo(T^) = ^ , (4) 
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with X ^ ^ Lagrange multipher fixed by the constraint on the average volume: 

V^{V)=J2Vp^iV) . (5) 

V 

EqH] is the - almost unavoidable - outcome of any statistical mechanics theory. The challenge is now to compute the 
number of microscopic states fl{V) associated with coarse-grained states which occupy a total volume V. To this end 
we must introduce some further assumptions as discussed in the next section. 



VII. SPACE PARTITIONS 



The phase-space volume fl{A) is the number of microscopic states which are classifiable under the same coarse 
grained state A. Specifically, ^(V) is associated with all the microscopic configurations $ = {ci, C2, c^} occupying 
a total volume V . Let us assume that all the cell-properties Ci are either completely determined by their volumes ut 
or they are independent from Vi . Let us also assume that these cells can have arbitrary volumes above the minimum 
value Vmim under the sole condition that the whole system must occupy a total volume V. In this case, the volume 
of the accessible phase-space can be computed exactly: 

1 /"^ /"^ /"^ (V - kij 

niV) = 1^^^ dv, j^^^^^ dv,.... 1^^^ dv,5{v, +v, + ... + v,-V)^ ^ A3.(r 1), ' (6) 

with A a constant analogous to the Debye length. Substituting into Eqdl and by using EqlS] we obtain X = ~ 
kvmin)/k and the probability 

(M\ f(M 1 \ ^ ~ ^min)^ ^ f 1 ^ ~ ^rain \ 

Poo{V) = f{V,k) = ———^ — - — — exp -fc-f^ , (7) 



r(fc) (y-y™„)'= "V y-Vr 

with Vmin — kvmin- The fuuctiou f{V,k) is the probability density function to find a packing of k elementary cells 
occupying a volume V, in a set of experiments where the average occupied volume is V. In a previous paper we 
reported a different derivation of f{V,k) from simple statistical arguments [s^l- Note that Eql6] is valid for any k 
and it holds even in the limit fc = 1. Indeed, the observable system can be any arbitrary sub-set of a larger system. 
Moreover, the experiment can be performed either on several different independent systems or - equivalently - on 
several non-interacting sub-sets of a large system. EqjT] is a Gamma distribution in the variable V — Vmin ] it is 
characterized by a 'shape' parameter k and a 'scale' parameter {V — Vmin) I k [33j . In the literature empirical fits 
with gamma distributions have been previously reported for two-dimensional Voronoi' networks [s^ ] , and for Voronoi 
partitions from random Poisson points in three dimensions [s^ . Interestingly the same distribution has been observed 
in statistical study of crumpled paper [36j . 

The average volume from the distribution /(V, k) coincides with V and the variance is 

2 {y ~ VminY 
" = 1 ' 

which is a useful relation to empirically evaluate k from a set of volume measurements: k — {V ~ Vmin)'^ /o'l- We will 
make use of this relation in Section HXl 



VIII. THE 'DEGREES' OF SPACE PARTITION 



Following Edward's ideas |4l.l37l|. a 'granular teniperature' can be inferred from an analogy with the thermodynamical 
relation P — l/{kBT) = d{Entropy)/d (Energy) [3,[33|, by susbstituting the volume to the role played by the energy 
in thermodynamical systems. In the formalism of this paper this can be written as: Pgr — X^^ = dS{Z)/dV. From 
Equations [2] and [5] we have: 



p-^=^=}L^ , (9) 

therefore the Edwards' compactivity x [li is the average free-volume per elementary cell, implying that the 
'granular temperature' is a measure of the kind and the degree of space-partition into elementary cells. For a given 
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system at a given packing fraction, the 'granular temperature' is univocally associated to the quantity k. By using 
Eqsl5]and[n] we have 

*^ * rmn 

The present theory provides therefore a very effective and practical way to compute Edward's compactivity from 
measures of volume fluctuations. 

In terms of the compactivity x the expression for the cell- volumes' distribution (Eq|7] ) can be written in the more 
compact form 

which recalls the equivalent expression in classical thermodynamics, with fl(V) given by Eql6] and ri(x) = 
Qqx'' exp{—V,nin/x)- The volume fluctuations within the ergodic set can be directly calculated from EqHand one 
can verify that the correct relation between compactivity and fluctuations is attained: 

x'^-{{V~{V)r) . (12) 
From this equation, substituting Eqs. [8]and[9l we obtain the following relation for the quantity k: 

k measures therefore the amount of volume that must be added to the system in order to increase of one 'granular 
degree' the compactivity. The analogous quantity for molecular gasses is: dEnergy /dTemperature, which is the 
specific heat. In analogy with ordinary thermodynamics this 'specific heat' is expected to be sensitive to changes in 
the system's internal properties. Ultimately, k counts the number of elementary cells, it is therefore a measure of the 
number of degrees of freedom. 

In this section we have shown that the present theory is consistent with all the equilibrium requirements in an 
ordinary statistical mechanics approach such as the molecular gas thermodynamics. Pushing forwards such analogy 
we can see that if the compactivity is the analogous of the temperature in a thermodynamics system, then Eql7]must 
correspond to the Maxwell-Boltzmann distribution for granular systems. 



IX. THEORY VALIDATION: VORONOI VOLUME DISTRIBUTION 



We have tested the validity of the present theory over a set of packings generated by means of a modified 
Lubachevsky-Stillinger algorithm [ssl. [33. |40| . The algorithm starts from random points in space and it makes them 
uniformly grow into non-overallping hard spheres with the sphere positions evolving in time according to Newtonian 
dynamics. The algorithm uses a cubic box and periodic boundary conditions. The simulation is ended when a jammed 
state with diverging collision rate is reached. Large expansion rates produce jammed configurations at low packing 
fractions whereas slower growth rates lead to larger packing fractions. For very slow rates, crystalline nuclei with large 
packing fractions are formed. In our simulations we used 10000 spheres with unitary diameters obtaining jammed 
configurations with packing fractions between 0.56 and 0.65 by varying the growth rate between 500 and 0.000001. 
We also generated non-jammed configurations in the range of packing fractions between 0.1 and 0.55 by keeping the 
growth rate at 0.001 and arresting the simulation once the desired packing fraction was reached. We also compare the 
theoretical predictions with experimental results concerning both acrylic beads in air and glass beads in water from 
the ANU database on disordered packings [4l| . Details on these experimental systems are reported in [2^, |3l, \^ . 

Within the framework of the present theoretical approach the packings generated by means of numerical simulations 
can be considered as prepared with a 'single trial'. The ergodic set for the jammed packings can be identified with the 
set of all configurations that can be achieved with a given growth rate. For the non-jammed packing configurations 
the ergodic set can be identified with the ensemble of non-overalpping spheres at a given packing fraction. In the case 
of dry acrylic and glass beads experiments the ergodic set is the set of all mechanically stable packings achievable with 
the specific preparation action [43|. Finally, we point out that, in the case of the glass beads in water, the preparation 
by means of a fluidized bed technique was specifically tuned with a sufficiently large number of repeated flows to 
ensure that the stationary state was reached jS^] . The fact that in all the above systems the statistical equilibrium is 
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reached is a 'working hypothesis' that will be supported, later in this section, by the very good agreement between 
the observations and the theoretical (equilibrium) predictions. 

Equation [7] is the main outcome of the present theory and it can be directly tested on the observations from 
experiments and computer simulations. The distribution f(V,k) predicted by Eq. [7] is valid for any aggregate of k 
elementary cells. Therefore, the same equation must describe both the volume fluctuations at the level of a single 
grain as well as the fluctuation at the level of the whole system. 

A. Local volume fluctuations 

Let us first study the system at local level investigating the volume fluctuations at the level of a single grain. For this 
purpose we must first implement a method to assign a portion of volume associated to each grain in the packing. A 
natural way to divide space into cells 'warped' around each grain is the Voronoi partition where cells are associated to 
the portion of space closest to a grain center respect to any other centre in the packing. Such a partition is particularly 
meaningful in the case of packing of equal spheres because in this case the Voronoi cell is always circumscribing the 
internal sphere. Extension to grains with more generic shapes and sizes can be implemented. 

Figure [1] shows the resulting distribution of the Voronoi volumes. One can see that such distributions span a very 
broad interval of volumes between 0.69 and 12 with large differences between different samples. Lines in FiglHare 
the plots of f{V,k) (EqlT]) with the parameter k calculated from k = {V — Vminf'/crt (Eql8|). One can see that all 
the distributions follow well the theoretical prediction. However, given the broad and different domains in which 
such distributions are defined, a better visualization of such agreement in necessary. This can be obtained by an 
appropriate shift and rescaling of the variable. Indeed, from Eq. [71 one can see that distributions characterized by 
similar values of k must result into similar behaviors when plotted as function of {V — Vmin)/{V — Vmin)- Figure 
shows the plot of all the distributions in function of such shifted-rescaled variable. We note that all the distributions 
from jammed configurations (bottom figure) collapse into very similar behaviors well described by /(V, k) with k 
in the range between 12.2 and 14.5. Un-jammed configurations (top figure) also reveal very good agreement with 
/(V, k) but in this case the parameter k spans a larger interval (9 < fc < 25). We also note that the distribution with 
crystalline inclusions (x symbols) deviates slightly from the predicted behavior of /(V, k) revealing the emergence of 
a second peak at small volumes (corresponding to the close-packed crystalline inclusions). 

The numerical results analyzed in the previous paragraph refer to rather 'ideal' packings made of non overlapping 
identical spheres. Remarkably, the sam e g ood agreement with f{V, k) (Eq[7]) was also found [32| in real experimental 
systems of acrylic spheres in air [2^ [s^, 113, ^^'^ glass spheres in water ^32] packed at different packing fractions 

by following different preparation protocols. Figure [3] shows that all the experimental distributions and the jammed 
configurations obtained numerically in the present work, collapse into very similar behaviors well described by /(V, k) 
with fc in a narrow range around k ~ 12. We note that, the probability distribution function is well reproduced in 
all the range of volumes; over eighteen different experimental systems; over eight different computer simulations; over 
a range of packing fractions between 0.55 to 0.64 and with a statistics performed over more than one million local 
configurations. 

The impressive fact of such an agreement is that these systems are very different (ideal Newtonian spheres, acrylic 
beads in air and glass beads in water) and they are prepared in very different ways (pouring, tapping, fluid flows, 
molecular dynamics shearing). The collapse of all these distributions around f{V,k ^ 12) suggests that there are 
universal properties that determine the packing conflgurations which are very little sensitive to the preparation 
method. Interestingly, we observe that un-jammed configurations are also very well described by f{V, k) but with 
parameters k which are different from the one of the jammed case. Figure |4] reports the estimates for the parameter 
k calculated from the relation k = {V — Vminf' 1(^1 (EqE]). The value at zero packing fraction (k = 5.586) was 
calculated analytically for random Poisson points in three dimensions [3^[4^. Figure [3] shows that, for non jammed 
configurations, the value of k increases almost linearly by increasing the packing fraction. Then the value of the 
parameter k collapses around 12 when the system gets into a jammed configuration. The inset in the figure shows 
that there are sizable differences between different systems and between the same system at different packing fractions. 
A careful look shows a rather neat variation in the value of k for the experimental systems around the packing fraction 
p ~ 0.6. Special properties associated to the global volume fluctuations at p ^ 0.6 was also observed in fluidized bed 
experiments by Schroter et al. Above the packing fraction ~ 0.645, the packings generate partially crystallized 

regions and the change in the kind of structural organization is reveled by a sharp drop in the value of k that eventually 
will tend to zero in the crystalline limit. We observe that the distribution /(y, fc) (Eq. [7]) with fc ~ 12 is also followed 
by the Voronoi volumes from simulations of granular packings reported in [46] . An equivalent collapse for a complete 
different system concerning for simulations of a polymer melt, water, and silica was reported by by Starr [47| where 
volume distributions were plotted vs. (v — {v))/av. It is straightforward to see that such a scaling coincides to the 
one proposed in this paper, indeed (w — {v))/a — Vk{{v — Vmin)/{.W) ~ Vmin) — !)• A careful look at the Starr's data 
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FIG. 1: (Color online), (top) Voronoi volume distributions from numerically generated equal-sphere packings at packing fractions in a 
range between 0.1 and 0.65. Different symbols refer to packings with different properties: jammed (o symbols), un-jammed (+ symbols) 
and packing with crystalline inclusions (x symbols). The lines arc plots of Eq[7]with the parameter k calculated from k = {V — Vmin)^/^v 
(EqlSjl. (bottom) The details of the jammed distributions at packing fractions between 0.55 and 0.65. 
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FIG. 2: (Color online). Voronoi volume distributions plotted vs. {V — Vmin) / {V — Vmin) ■ (top) All distributions: jammed (o), un-jammed 
(+) and packing with crystalline inclusions (x). (bottom) The details of the jammed data show that they all collapse onto distributions 
in good agreement with Eq[7]with very similar values of k (between 12.2 and 14.5). 
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FIG. 3: (Color online). Comparison between the experimental Voronoi volumes distributions for all the 18 experimental samples with 
acrylic and glass beads in air and water |32|| and all the numerical jammed configurations (o) obtained by using the modified Lubachevsky- 
StiUinger algorithm [s^lsl, [Z^l. Experiments with dry acrylic beasd in air have symbols: A (<), B (□), C (*), D (o), E (A), F (v)- 
Experiments with glass beads in water have symbols +. When plotted vs. {V — Vmin)/iV — Vmin) a^U the data collapse to a common 
behavior which is described well by f{V,k = 12) (black line). 

indicates a good agreement with /(V, k ^ 20). 

B. Global volume fluctuations of the whole system 

We have seen in the previous discussion that the volume fluctuations at the level of a single grain are very well 
described by f{V, k) and we have observed that the value of the parameter fc is a good indication of the system state 
(un-jammed, jammed, poly-crystalline). Let us now verify if the present theory can also correctly predict the volume 
fluctuations at the level of the whole sample. This requires the analysis of the distribution of the volumes of large 
number of equivalent samples made by means of the same procedure. Differently from the single-grain case here we 
can only study jammed configurations. We have prepared such samples by using the modified Lubachevsky-Stillinger 
algorithm [SSl [39l, |40| generating a large number of jammed packings containing different numbers of grains (between 
N — 200 to TV = 10000) and prepared at different growth rates (between 0.0005 and 1000). For each growth rate and 
size we run at least 200 independent simulations and we measure the total volume occupied by the jammed spheres. 
Three examples of the resulting statistical distributions of volumes are shown in Fig[51 Again we find that Eq[7|predicts 
extremely well the observed fluctuations. The agreement is obtained without using any adjustable parameter. Indeed, 
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FIG. 4: (Color online). Behavior of k calculated from k = (V — Vmin)^ I (Eq[8ll vs. packing fraction p. This plot refers to the statistics 
on local configurations computed over the fluctuations of the Voronoi cells volumes. Neat changes in k are observed when spheres first 
jam (p ~ 0.55) and when crystalline nuclei start to form (p > 0.65). The '+' refers to simulations of un-jammed packings of spheres; the 
'o' are jammed packings simulated by using different growth rates and 'x is a jammed sample with crystalline inclusions; the symbols '!>' 
correspond to 12 experiments with glass spheres in water prepared by means of fluidized beds technique 32j ; finally, the 6 symbols 
□, it, o, A , V ' refer to 6 experiments with acrylic beads in air prepared by means of different methods )28i . i43 . 14^ . |44| . 



to calculate the parameter fc we compute the average volume ly) = V , the variance cr^ — 

_ (y)2 and we use 

k = [y — Vrnin)'^ /<^v (EqlS)). In this case, Knm is fixed to the minimum volume attainable by a packing of N equal 
spheres with unit diameter, which is the volume occupied by a fee crystalline configuration: i.e. N/ \/2 (for spheres 
with unitary diameters). 

One can note that the distribution tends to narrowing during compaction. The quantity k is changing with the 
packing fraction and, for example, in the three examples in FigOwe measure: k = 2.6N at p — 0.564, k — AAN 
at p = 0.616 and k — 8.8N a.t p — 0.637. More generally we compute the average volume and the variance for a 
very large set of data over a range of packing fractions within 0.55 and 0.66. The resulting values of the parameter 
k rescaled by the number of spheres in each sample N are plotted in FiglH We note that the data from samples 
with various sizes (between = 200 to = 10000) all fall on similar values when rescaled by N indicating that the 
parameter k scales linearly with the number of spheres. The rescaled parameter k/N takes values a little below 3 at 
the loose packing limit {p ~ 0.55). Then the value of the parameter grows consistently when the packing densifies 
and it reaches a peak at the random close packing limit p ~ 0.64. Afterwards, k/N drops sharply and eventually goes 
to zero when the packing becomes crystalline. The small samples with 200 spheres show a slightly different behavior 
with the peak occurring before p ~ 0.64. We verified that this is a finite-size effect due to the formation of some 
poly-crystalline samples. 

It is therefore clear that the quantity k is very sensitive to the changes in the system's internal structure and the peak 
at the random close packing limit indicates that the packing must reorganize in order to proceed with compaction. A 
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FIG. 5: (Color online). Fluctuations of the whole sample volume (symbols) from large number of repeated simulations of jammed 
packings. The distributions are from samples each containing n = 1000 spheres and generated with growth rates 0.01 (p = 0.637), 0.05 
(p = 0.616) and 1000 (p = 0.564). The lines are the theoretical prediction from Eq[7] They use no adjustable parameters: the parameter 
k is calculated from k = {V — Vminf' I (EqO and Vmin = 



similar discontinuity at both the loose and closed packing limits was observed in the study of k from the local Voronoi 
volumes statistics (FigU]). However, the values of the local quantity (fc) and the rescaled global quantity (k/N) are 
different. This should not be of any surprise. Indeed, k counts the number of elementary cells contributing to the 
system's volume. But, some elementary cells are sheared among neighboring Voronoi cells, and therefore the number 
that contributes to the volume of N Voronoi cells is not necessary N times the number that contributes to one single 
cell. Furthermore, it is known that volume-correlations play an important and unusual role in these systems |,4S |. It 
has been shown both in computer simulations and experiments [Hj^^lsOl that the scaling of the volume fluctuations 
is not a simple linear function of the portion of sample investigated and there are evidences of a dependence on the 
preparation procedure and on the packing fraction. Although not surprising, such different values and behaviors of 
the local and global k have important implications. For instance this implies that the granular temperature measured 
at local level is different from the one measured from the global fluctuations. In this respect the 'proper' granular 
temperature is the one associated to the global volume fluctuations which we have shown is an intensive quantity since 
k is scaling linearly with TV (at least for samples with more than 200 spheres, FiglS]). The 'temperature' associated 
with the local measure of k is an 'effective granular temperature'. The fact that they predict so accurately the volume 
fluctuations both at local and at global level indicates that the effect of correlations and local geometrical constraints 
can be treated within a statistical mechanics framework by simply tuning the value of k. 
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FIG. 6: (Color online). The parameter k divided by the number of spheres in the sample (A'^) is plotted versus the average packing 
fraction p = Tv/iGV) for several samples with total number of spheres between 200 and 10000. The value of k is computed from k/N = 
1/N(V — V,„i,i)/cr^ ~ ^ "^^"^ I'^'i ^ la^rge set of samples prepared by repeating a large number of simulations using the same growth 
rate for each average packing fraction. The full line is a guide for eyes showing the sharp transition occurring around 0.645 for samples 
with N > 700, the dotted line highlights that small samples (A'^ = 200) start the decay at smaller packing fractions around p ~ 0.635. 



X. CONCLUSIONS 

In this paper we have shown that a statistical mechanics theory for granular matter can be always deductively 
constructed from two simple physical assumptions: 1) the microscopic system's structure can be completely described 
in terms of a finite set of quantities; 2) after performing a given experiment, the probability to find the system in a 
certain static state depends only on the previous state (and on the experiment itself). These two simple assumptions 
lead to a complete statistical mechanics description for these systems where both transient and equilibrium properties 
can be treated. We show that an expression for the volume fluctuations can be explicitly obtained (Eq|3]) by introducing 
the following three additional assumptions: i) the microscopic state can be encoded in terms of the properties of a 
set of local cells with volumes Vi\ ii) the cells properties are either fully described by their volumes Vi or they are 
independent from Vi\ iii) any space partition with cells of volumes vi > Vmin satisfying the condition '^^Vi = V 
corresponds to an attainable packing. A comparison with a large number of experiments and computer simulations 
shows remarkably good agreements with the theoretical predictions in a wide range of packing fractions and for several 
differently-prepared systems (Section HXl and Figs. [TJ O and O . 

We have identified a quantity, fc, which is very sensitive to the system's structural organization. It results that such 
quantity is the equivalent of the speciflc heat in ordinary thermodynamics. We have demonstrated that k can be easily 
computed empirically (Eq|5]) by measuring the average volume and the variance in a set of trials. A large number of 
measures from both the local statistics at the level of the single grain and at the global level of the whole sample show 
the signatures of clear structural transitions marked by peaks and discontinuities in the values of k occurring at both 
the random loose packing and at the random close packing limits (Figs |4] and [6]). Similar peaks in the specific heat 
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are observed at phase transitions in thermodynamic systems. Understanding whether proper thermodynamical phase 
transitions are occurring at the he random loose packing Umit and the random close packing limit is the subject of 
current investigations. 

We have demonstrated that the present theory can predict both the volume fluctuation at the level of a single grain 
(Voronoi' cells) and at the level of the whole system (see Figs. [3] and [5]). The remarkable quantitative agreement 
between theory, experiments and simulations at both local and global levels is very encouraging demonstrating that 
statistical mechanics is a powerful tool to study and predict the properties of complex materials. We have discussed 
that the local effects of correlations and geometrical constraints can be taken into account by using an 'effective' 
local granular temperature. Further work must be devoted to better understand the link between local and global 
properties and identify the length-scale above which the system becomes extensive. An important test to establish the 
correct granular temperature might use an independent measure from the fluctuation-dissipation theorem as discussed 
in 0, The use of non extensive forms of entropy instead of the classical Gibb's entropy might be right approach to 
formulate a theory which consistently describes both local and global levels in terms of a unified generalized granular 
temperature. 
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